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ABSTRACT 

During a lake-wide survey in 1963, measurements were made of depth, % 
organic matter, mixed sediment oxygen uptake rate, and sediment core oxygen 
uptake rate. 

Depth had a poor, negative correlation to other variables except % organic 
matter, which correlated positively. Percent organic matter averaged 6.396%, 
and contours followed current patterns. It also had a strong correlation with 
mixed sediment oxygen uptake rate, which averaged 318.6 ]jg 02/gm sediment/5 
minutes. Core oxygen uptake rate, averaging 4.987 yg 02/cm^/hour, correlated 
strongly between incubation durations and moderately well with mixed sediment 
oxygen uptake rate. Cluster analysis considering all variables had an 
acceptable cophenetic correlation of 0.7247. The ratio of mixed sediment 
oxygen uptake rate to % organic matter indicated the oxidizability of sediments 
and was higher toward the west. Clusters and ratio values were used to define 
lake regions of similar sediment characteristics. Regression models of mixed 
sediment oxygen uptake rate versus % organic matter were computed for each 
cluster analysis region and ratio region. Slopes for linear regression models 
generally increased toward the west, indicating the increasing oxidizability of 
the sediments . 
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INTRODUCTION 

The eutrophication of Lake Erie is indicated by a hypo limnetic 
dissolved-oxygen depletion during summer stratification. This has been of 
major interest ever since the possibility and eventual reality of anoxic 
conditions in the central basin were noted in the late 1950s. 

The history of this problem is well documented (Dobson and Gilbertson 
1972, Carr 1962, Carr et al. 1965), despite the absence of synoptic data prior 
to the 1960s. Significant changes in benthic species composition and abundance 
showed the consequences before anoxic conditions were directly observed. Since 
first noted in 1929, this oxygen depletion has become more severe, affecting a 
larger area more quickly after stratification (Carr 1962, Dobson and Gilbertson 
1972). 

Anoxic conditions in the hypolimnion present several problems, as outlined 
by Burns and Ross (1972b). The adverse effects upon benthic fauna are obvious, 
as seen in the elimination of the mayfly (H exagenia ) population, as well as in 
a shift toward more tolerant oligochaetes. A possible consequence, and of 
greater economic importance, is the replacement of more valuable fish such as 
blue pike and Cisco by yellow perch and smelt as the basis for the lake's 
fishery. As the anoxic region became more widespread, difficulties were 
incurred by municipalities which use the lake as a water source due to the 
objectional taste and smell of anoxic water. Also of major concern is the 
increased dissolution of nutrients from sediments in the absence of oxygen, 
creating the possibility of a self-fertilization cycle. 

Lake Erie is divided into three distinct basins that behave differently, 
primarily due to differences in bathymetry (Figure 1). 

The western basin is extremely shallow, receives a large influx of 
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nutrients and water via the Detroit, Maumee, and Raisin Rivers, and rarely 
stratifies (Schelske and Roth 1973). Although the western basin is highly 
productive with a large amount of organic matter in the sediments, the mixing 
due to wind and river inflow is usually sufficient to prevent anoxic conditions 
(Potos 1970). 

Considering its vast area of 16,300 km,^, the central basin is relatively 
shallow, with an average depth of 18.3 meters. In spite of the long wind 
fetch, the basin will stratify with a thermocline occurring near 17 meters. 
This results in an average hypolimnion thickness of only 2.5 meters (Burns and 
Ross 1972b, Dobson and Gilbertson 1972). After isolation from the atmosphere, 
the hypolimnion is probably depleted of oxygen primarily by the sediments, due 
to the large sediment surface area to hypolimnion volume ratio (Lucas and 
Thomas 1971) . 

The deep eastern basin does not suffer from anoxia, even though it acts as 
a settling area for organic matter, due to the large hypolimnion volume and the 
low oxygen demand exerted by the sediments . 

The longstanding belief has been that the oxygen demand in hypolimnetic 
water is due primarily to decomposition and bacterial oxidation within the 
water column of sinking organics , predominantly dead phytoplankton (Kleveno 
et al. 1971). However, the BOD of the overlying water is seen to be 
insufficient to cause anoxia in the length of time observed. The initial 
purpose of this study and other recent research was to find alternative 
mechanisms of oxygen depletion to explain this discrepancy. An attempt has 
also been made to quantify depletion rates. 

Potos (1970) attributed oxygen depletion to respirational requirements of 
bacteria in the degradation of organic matter, primarily sedimented plankton, 
and to demand by reduced chemicals when exposed to overlying water. Citing the 



occurrence of unusually high hypo limnetic current velocities of 61 cm/sec, 
Potos suggested that sediments can thereby be resuspended and a chemical demand 
created. 

This theory was supported by Project Hypo, an intensive study of the 
central basin in 1970 (Burns and Ross 1972a). Blanton and Winklhofer (1972) 
reported average hypolimnetic current velocities between 2.2 and 5.6 cm/sec, 
but noted velocities of up to 98.5 cm/sec. They also provided additional 
information on the hypolimnion' s dynamic nature, a detailed physical 
description, and documentation of volume changes. Sediment oxygen demand 
measured at 0.31 gm 02/m^/day could account for the observed depletion rates 
but was not assumed to be the only parameter involved. 

A mechanism postulated by Kleveno et al. (1971) for sediment oxygen 
depletion involving benthic algae was tested in Project Hypo. The results 
(Braidech et al. 1972) showed that two species of algae, Tribonema utriculon 
and Qedogonium sp., remained viable after settling, but their contribution of 
oxygen was masked by stronger oxygen-demanding physical and chemical phenomena. 

The conclusion that benthic algae were not of major consequence was also 
reached by Lucas and Thomas (1971). They did, however, observe a diurnal 
oxygen fluctuation. In situ sediment oxygen demand was calculated as 0.40 gm 
02/m^/day. 

Burns and Ross (1972c) determined an aerial uptake rate of 0.39 gm 
02/in^/day. This value is considered most representative of Project Hypo 
results since it involved a basin-wide survey. 

Burns (1976) asserted that oxygen depletion rates could serve as important 
indicators of the trophic state of Lake Erie if rates could be measured to an 
accuracy of 3% over a time scale of 5 years. Burns modeled hypolimnetic 
reoxygenation through oxygen budgets and calculated an aerial depletion rate of 



0.43 gm 02/ni /day and a voliimetric depletion rate of 0.13 gm 02/ni'^/day (3.9 mg 
C)2/L/inonth) . 

Dobson and Gilbertson (1972) used historic changes in oxygen concentration 
through the summer seasons to calculate oxygen depletion rates. They clearly 
show a progressive increase in the oxygen depletion rate over the last 40 
years. Based on a length of stratification of 110 days, a volumetric depletion 
rate of 3.0 mg 02/L/month will result in anoxic conditions. This critical rate 
was reached about 1960, and Carr (1962) reported anoxic conditions over 3600 
km^ in 1960. 

This study attempts to elicit the relationship between organic matter and 
sediment oxygen uptake rate as well as to perceive the relative importance of 
this in overall hypolimnetic depletion. 

PROCEDURKS 

After a preliminary study in the summer of 1962 to assess sampling 
methods, a lake-wide survey was completed in 1963 by the Bureau of Commercial 
Fisheries. Employing the research vessels Cisco and Kaho , cruises were made of 
the western and central basins in July and of the eastern basin in September. 
At each station, (Table 1, Figure 2) depth was measured and an Ekman grab 
sample taken. Additionally, certain physical and chemical properties of the 
overlying water were measured at some stations. 

Imtaediately after the Ekman grab was retrieved, about 0.25 gm (dry weight) 
of the top few millimeters of sediment was skimmed off, placed in a 32.5-mL 
D.O. bottle with water of known saturated oxygen concentration, stoppered by an 
oxygen electrode, and kept in suspension using a magnetic stirrer. Oxygen 
consumption for the mixed sediment was monitored over 5 minutes. These values 



TABLE 1. Sampling locations and dates. 





STATION # 


LATITUDE 


LONGITUDE 


DATE 


1. 


41035' 


45" 


82050 " 


00" 


11 9/63 


2. 


41039 ' 


15" 


82052 " 


30" 


11 9/63 


3. 


41043" 


55" 


82053 " 


00" 


11 9/63 


4. 


41048" 


00" 


82053 " 


55" 


11 9/63 


5. 


41052 ' 


45" 


82049 " 


10" 


11 9/63 


6. 


41057' 


05" 


82044 " 


10" 


11 9/63 


7. 


41052 ' 


45" 


82044 " 


40" 


11 9/63 


8. 


41048' 


30" 


82045" 


10" 


11 9/63 


9. 


41043" 


00" 


82 046" 


00" 


11 'i/di 


10. 


41038" 


25" 


82045 ' 


50" 


7/10/63 


11. 


41033" 


40" 


82045 " 


45" 


7/10/63 


12. 


41033" 


45" 


82042 


30" 


7/10/63 


13. 


41032 " 


40" 


82 040 " 


00" 


7/10/63 


14. 


41036" 


00" 


82038' 


35" 


7/10/63 
9/ 5/63 


15. 


41040 " 


15" 


82035 


20" 


7/10/63 


16. 


41044" 


45" 


82034 


45" 


7/10/63 


17. 


41048" 


25" 


82034 


40" 


7/10/63 


18. 


41053" 


00" 


82034 ' 


35" 


7/10/63 


19. 


41057" 


45" 


82034 


30" 


7/10/63 


20. 


42O00" 


20" 


82026 


50" 


7/11/63 


21. 


41056" 


00" 


82026 


25" 


7/11/63 


22. 


41051 " 


40" 


82026 


05" 


7/11/63 


23. 


41047 " 


00" 


82026 


45" 


7/11/63 


24. 


41042' 


50" 


82 02 7 


50" 


7/11/63 


25. 


41038" 


25" 


82029 


00" 


7/11/63 


26. 


41034" 


00" 


82 030 


'10" 


7/11/63 


27. 


41029 


'50" 


82031 


'10" 


7/11/63 


28. 


41028' 


00" 


82022 


'15" 


7/11/63 


29. 


41026 


00" 


82027 


'45" 


7/11/63 


30. 


41025 


30" 


82032 


'20" 


7/11/63 


31. 


41032 


05" 


82022 


'30" 


7/12/63 


32. 


41036 '20" 


82022 '30" 


7/22/63 



(continued) 



TABLE 1. (continued) 





STATION # 


LATITUDE 


LONGITUDE 


DATE 


33. 


41041' 


30" 


82 02 2 ' 


30" 


7/12/63 


34. 


41045" 


00" 


82 021 ' 


40" 


7/12/63 


35. 


41049" 


15" 


82O20' 


25" 


7/12/63 


36. 


41053' 


20" 


82 01 9 " 


20" 


7/12/63 


37. 


41057' 


50" 


82 01 8' 


10" 


7/12/63 


38. 


42O02' 


00" 


82017 ' 


05" 


7/12/63 


39. 


42O05' 


55" 


82 01 4' 


15" 


7/12/63 


40. 


42O08' 


10" 


82O08' 


50" 


7/12/63 


41. 


42O10 


45" 


82 00 4' 


10" 


7/12/63 


42. 


42O09 


30" 


81056' 


30" 


7/12/63 


43. 


42O05 


20" 


8I058 


10" 


7/13/63 


44. 


42O01 


00" 


81 059 


40" 


7/13/63 


45. 


41057 


00" 


82 001 


15" 


7/13/63 


46. 


41052 


55" 


82 002 ' 


40" 


7/13/63 


47. 


41048 


30" 


82 00 4 


10" 


7/13/63 


48. 


41044 


30" 


82 005 


45" 


7/13/63 


49. 


41040 


'20" 


82 00 7 


'15" 


7/13/63 
9/ 5/63 


50. 


41036 


'00" 


82O08 


40" 


7/22/63 


51. 


41032 


'15" 


82 00 3 


'30" 


7/13/63 


52. 


41031 


'50" 


82 010 


'15" 


7/13/63 


53. 


41035 


'15" 


81 05 7 


'15" 


7/15/63 


54. 


41039 


25" 


81 055 


40" 


7/22/63 


55. 


41043 


'40" 


81 054 


25" 


7/15/63 
9/ 5/63 


56. 


41047 


'35" 


81 052 


'40" 


7/15/63 


57. 


41051 


'50" 


81 051 


'00" 


7/15/63 


58. 


41056 


'00" 


81 049 


'30" 


7/15/63 


59. 


42 000 


'00" 


81 04 7 


'55" 


7/15/63 


60. 


42O04 


'10" 


81046 


'15" 


7/15/63 


61. 


42O08 


'45" 


81 044 


'40" 


7/15/63 


62. 


42011 


'25" 


81 049 


'00" 


7/15/63 


63. 


42013'15" 


81 044 '00" 


7/16/63 



(continued) 



TABLE !• (continued) 





STATION # 


LATITUDE 


LONGITUDE 


DATE 


64. 


42Ol2'00" 


81034.55" 


7/16/63 


65. 


42O10'40" 


81O26'20" 


7/16/63 


66. 


42O04'20" 


81O29'20" 


7/16/63 


67. 


4IO58'10" 


81O32'10" 


7/16/63 


68. 


41052'00" 


81O35'10" 


7/16/63 


69. 


41045'50" 


81O38'10" 


7/16/63 
9/ 5/63 


70. 


41039 '40" 


81041'10" 


7/17/63 


71. 


41033 '30" 


81044'10" 


7/17/63 


72. 


41035'50" 


81 038 '10" 


7/17/63 


73. 


41041'30" 


81O33'40" 


7/17/63 


74. 


41047'00" 


81028*50" 


7/17/63 
9/ 5/63 


75. 


41052'30" 


81O24'20" 


7/17/63 


76. 


4lO58'00" 


81019'30" 


7/17/63 


77. 


42O02'00" 


81016'10" 


7/17/63 


78. 


42O06'45" 


81012'30" 


7/17/63 


79. 


42O03'15" 


81004'40" 


7/17/63 


80. 


41059 '40" 


80O57'30" 


7/17/63 


81. 


41056'15" 


80O50'20" 


7/17/63 


82. 


41054 '00" 


8lO00'00" 


7/17/63 
9/ 5/63 


83. 


41057 '45" 


81 000 '00" 


7/18/63 
7/21/63 


84. 


42O04'10" 


8lO00'00" 


7/18/63 


85. 


42Ol0'25" 


8lO00'00" 


7/18/63 


86. 


42O17'00" 


81 000 '00" 


7/18/63 


87. 


42O23'10" 


SIOOO'OO" 


7/18/63 


88. 


42O29'10" 


81 002 '55" 


7/18/63 


89. 


42034 '20" 


81O07'30" 


7/18/63 


90. 


42033 '45" 


81016'00" 


7/18/63 


91. 


42027'45" 


81018'35" 


7/19/63 


92. 


42 021 '50" 


81O21'20" 


7/19/63 



(continued) 



TABLE 1. (continued) 





STATION # 


LATITUDE 


LONGITUDE 


DATE 


93. 


42015' 


50" 


81 02 3 


55" 


7/19/63 


94. 


42015" 


45" 


81 031 


10" 


7/19/63 


95. 


42022' 


00" 


81 03 8 


'45" 


7/19/63 


96. 


42026' 


15" 


81 032' 


30" 


7/19/63 


97. 


42O30 


40" 


81 02 6 


'00" 


7/19/63 


98. 


42 035 


00" 


81 020 


00" 


7/19/63 


99. 


42034 


55" 


80O56 


'30" 


7/20/63 


100. 


42032' 


40" 


80 048 


40" 


7/20/63 


101. 


42026 


20" 


8O046 


15" 


7/20/63 


102. 


42O20 


25" 


80 043 


'20" 


7/20/63 


103. 


42018 


00" 


8O032 


'40" 


7/20/63 


104. 


42014 


15" 


8O04O 


'10" 


7/20/63 


105. 


42 oil 


45" 


8O033 


00" 


7/20/63 


106. 


42O08 


15" 


80 037 


'20" 


7/20/63 


107. 


42O05 


15" 


8O033 


15" 


7/20/63 


108. 


42 002 


50" 


80 030 


'20" 


7/20/63 
9/ 6/63 


109. 


42 002 


'15" 


8O034 


'50" 


7/21/63 


110. 


42O02 


'40" 


8O042 


'00" 


7/21/63 


111. 


42O06 


10" 


80 049 


'45" 


7/21/63 


112. 


41055 


45" 


81 010 


•40" 


7/21/63 


113. 


41051 


'20" 


81017 


'40" 


7/21/63 


114. 


41034 


'50" 


81 051 


'20" 


7/22/63 


115. 


42O09 


'00" 


80O27 


'50" 


9/ 6/63 


116. 


42 01 5 


'20" 


80O25 


'15" 


9/ 6/63 


117. 


42019 


50" 


80O23 


'20" 


9/ 6/63 


C 9. 


42029 


'45" 


8OOI8 


'45" 


9/ 6/63 


118. 


42023 


'00" 


80 01 8 


'45" 


9/ 6/63 


CIO. 


42016 


30" 


8OOI8 


'40" 


9/ 6/63 


119. 


42O10 


'00" 


80 01 8 


'45" 


9/ 6/63 


120. 


42O10 


'30" 


8OOII 


'35" 


9/ 7/63 


121. 


42013 


'40" 


80O15 


'40" 


9/ 7/63 


122. 


42019 


'25" 


8O012 


'50" 


9/ 7/63 


123. 


42025 


00" 


80O09 


'50" 


9/ 7/63 


124. 


42031 '00" 


80O07'00" 


9/ 7/63 



(continued) 



TABLE 1. (continued) 





STATION # 


LATITUDE 


LONGITUDE 


DATE 


C25. 


42021" 


00" 


8OOO6" 


20" 


9/ 7/63 


C 4. 


42015" 


00" 


8OOO8" 


00" 


9/ 7/63 


125. 


42^12 ' 


00" 


80O05 


00" 


9/ 7/63 


C 5. 


42017 " 


15" 


79058" 


30" 


9/ 8/63 


C 7. 


42O30' 


00" 


79055" 


50" 


9/ 8/63 


126. 


42 032 ' 


45" 


80 002 " 


25" 


9/ 8/63 


127. 


42034" 


25" 


80 001 


30" 


9/ 8/63 


128. 


42039" 


00" 


80O13" 


00" 


9/ 8/63 


129. 


42043" 


00" 


80O13 


00" 


9/ 8/63 


130. 


42 046" 


15" 


80O12 


30" 


9/ 8/63 


131. 


42042 ' 


55" 


80O07 


00" 


9/ 9/63 


132. 


42039" 


00" 


80O02 


00" 


9/ 9/63 


133. 


42034" 


30" 


79056 


05" 


9/ 9/63 


134. 


42030 


00" 


79050 


00" 


9/ 9/63 


135. 


42023' 


00" 


79050 


00" 


9/ 9/63 


136. 


42016 


10" 


79050 


00" 


9/ 9/63 


137. 


42025 


00" 


79038 


00" 


9/ 9/63 


138. 


42030 


00" 


79024 


00" 


9/ 9/63 


139. 


42034 


20" 


79023 


'45" 


9/10/63 


140. 


42034 


'20" 


79033 


20" 


9/10/63 


141. 


42034 


20" 


79042 


'30" 


9/10/63 


142. 


42038 


45" 


79050 


'00" 


9/10/63 


143. 


42043 


00" 


79057 


'30" 


9/10/63 


144. 


42043 


00" 


79047 


'00" 


9/10/53 


145. 


42043 


'00" 


79038 


'30" 


9/10/63 


146. 


42048 


'45" 


79036 


'00" 


9/10/63 


147. 


42048 


'15" 


79028 


'10" 


9/11/63 


148. 


42043 


'20" 


790I8 


'30" 


9/11/63 


149. 


42048 


'25" 


79009 


'00" 


9/11/63 


150. 


42049 


'45" 


79000 


'20" 


9/11/63 


151. 


42046 


'15" 


78057 


'00" 


9/11/63 


152. 


42041 


'20" 


790O6 


'35" 


9/11/63 


153. 


42035 


'00" 


79015 


•00" 


9/11/63 


154. 


42038 


'00" 


79023 


•40" 


9/11/63 


155. 


42O38'00" 


79032*00" 


9/11/63 
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are referred to as mixed sediment oxygen uptake rates. The sediment was then 
filtered from the water, dried at 60^C, and weighed. Organic content was 
determined by ashing at 600^C. 

A separate portion of the surficial sediment was centrifuged to the bottom 
of several 50-mL vials; the vial sides were washed and recentrifuged, and the 
supernatant decanted. A known volume and oxygen concentration of oxygen- 
saturated water was added and the samples incubated for 2, 4, and 6 hours at 
20^C. Incubation water was then siphoned with a syringe into 5-mL glass- 
stoppered vials used as micro D.O. bottles. Oxygen consumption values 
determined by this procedure are referred to as core oxygen uptake rates. 
Possible oxygen contamination was not considered to be as critical as for other 
procedures, since this would only further minimize an intentionally 
conservative measure of oxygen uptake. 

The Winkler method was used for all D.O. determinations, with the micro 
D.O. bottles having 2 mL titrated with 0.0025 N sodium thiosulfate. 

Analysis using the University of Michigan's MIDAS statistics program was 
done on only those stations with complete data for depth, % organic matter, 
mixed sediment oxygen uptake, and core oxygen uptakes (Table 2). Variables 
were plotted by station location and contours drawn. Correlation coefficients 
were calculated between all variables (Table 3) and, based on these results, 
mixed sediment oxygen uptake rates were plotted against % organic matter. 

In order to further stratify the data, several cluster analyses were done 
using all variables. Additionally, the ratio of mixed sediment oxygen uptake 
to % organic content was calculated for each station. Regions determined by 
clustering and by similarity of ratio values were then used to plot % organic 
matter versus mixed uptake rate separately. Least square regression models 
were calculated for each region. 
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Scatter plots were also obtained by region for depth and % organic matter 
against core oxygen uptake and for depth versus mixed sediment oxygen uptake. 

RESULTS 

Sampling depths ranged from 7 to 62 meters (Table 4), with an average of 
20.5 meters, which is representative of the major portion of the central basin 
(Figure 1). Depth had a poor negative correlation to all of the other 
variables except % organic matter which had a positive correlation (Table 3). 
Scatter plots of depth versus core and mixed sediment oxygen uptake rates also 
indicated a lack of strong correlation. 

Percentage of organic matter varied from 0.25 to 10.76% and averaged 
6.396% (Table 4). This reflects a general trend away from shore (Figure 3), 
with the greatest values in the deeper basins and a scouring effect lowering 
values near shore, as well as off Pt. Pelee and Pte. Aux Pins. There is, for 
the most part, a paucity of organic content in sediments of the eastern third 
of the central basin. The south shore of the central basin is characterized by 
higher values offshore, extending from Sandusky, Ohio, to east of Erie, 
Pennsylvania. Average values are seen in the western basin, while the eastern 
basin exhibits higher values in deep regions and very low values in the far 
eastern end. Percentage of organic matter showed good correlation (0.7213) 
with mixed sediment oxygen uptake and less correlation to core oxygen uptake 
rates and depth (Table 3). 

Mixed sediment oxygen uptake had an equally broad range, from 8. to 754. 
]sg 02/gin sediment/5 minutes, averaging 318.6 yg 02/gm sediment/5 minutes (Table 
4). Contours (Figure 4) followed those of % organic matter fairly closely, 
with the highest uptake occurring in regions of high organic content. For a 
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given organic content, however, uptake was highest toward the west. This 
relationship can be seen by comparing mixed sediment oxygen uptake rates of 
400, 500, and over 600 ]ig 02/gia sediment/5 minutes for the eastern, central, 
and western basins, respectively, all in areas of about 9% organic matter. 
Compared with % organic matter, mixed sediment oxygen uptake rates had less 
correlation to the other variables (Table 3). 

Core oxygen uptake rates ranged from 0.36 to 16.19 yg 02/cm^/hr for 2-hr 
incubation but became more conservative with longer incubation durations . 
Averaged core oxygen uptake rates for each station ranged from 0.57 to 8.33 yg 
02/cm2/hr, with a mean of 4.987 \xg 02/cm^/hr (Table 4). Core oxygen uptake 
rates correlated only moderately well (0.4684-0.5471) with mixed sediment 
oxygen uptake rates, and poorly with other variables (Table 3). There was, 
however, excellent correlation (0.7366-0.9429) between uptake rates of various 
incubation times. Contours of average core oxygen uptake rates (Figure 5) were 
fairly consistent with contours of % organic matter and mixed sediment oxygen 
uptake rate, but scatter plots of core oxygen uptake rates versus % organic 
matter and depth showed no relationship between them. 

Cluster analysis is a statistical technique done by computer in which 

similarity of stations is determined by comparing relative proximity after 

plotting stations in N-dimensional space, where N is the number of variables 

characterizing each station. The cluster analysis that yielded the best 

grouping, or highest cophenetic correlation (0.7247), was standard Euclidean 

with averaged distances. While this optimal case exhibits only fair 

clustering, the groupings thus determined did fall into distinct regions 

(Figure 6) within the lake, with a few excluded stations that did not exhibit 

similarity to any others. Cluster regions had similar sediment characteristics 

based on depth, percent organic matter, mixed sediment oxygen uptake rate, core 
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oxygen uptake rates, and ratio of mixed sediment oxygen uptake rate to percent 
organic matter. 

Ratio values, a variable dependent on % organic and mixed uptake rate, 
rarely had a value greater than one. This value was computed by dividing the 
mixed sediment oxygen uptake rate by 100 times the % organic matter for each 
station. These values, therefore, show the degree of oxidizability of 
sediments or, in other words, the oxygen demand of sediment with a certain % 
organic matter. Curiously, ratio values (Figure 7) showed a low correlation 
(0.2646) to % organic matter (Table 3). These values were used to divide the 
lake into regions of similar sediment oxidizability (Figure 8). Highest 
values, >0.70, were found in the western basin and the southwest corner of the 
central basin (Regions 1 and 2). These graded into those between 0.70 and 0.40 
in Region 3, while values in the high 0.60s were found in Region 4. Regions 5 
and 6 were both characterized by ratio values generally less than 0.40. 

Scatter plots of mixed sediment oxygen uptake versus % organic matter 
based on grouping by cluster analysis and ratio similarity exhibited either a 
tight clumping or a significant regression. Those regions where the regression 
was not meaningful, significance >0.05, do not have the least squares model 
indicated. Regressions for cluster analysis regions (Table 5) had lower R 
values (0.46873-0.77617) and higher significance indices than those for ratio 
analysis regions (Table 6) (0.70091-0.97991), indicating less significant 
cluster analysis regressions. Slopes of the linear regression models decreased 
successively (95.71, 35.17, 27.18, 23.86) through Regions C, A, D, and B, 
respectively, for the cluster analysis, and (125.77, 79.97, 68.75, 57.70, 
47.56, 41.06) through Regions 1, 2, 4, 3, 6, and 5, respectively, for the ratio 
analysis . 
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TABLE 5. Oxygen uptake of mixed sediment related to percent organic matter for 
regions determined by cluster analysis in Lake Erie 1963. Y = yg O2 consiamed/ 
gm sediment/5 minutes; X = % organic matter. 



Regions A and Al. Stations 2-6, 8-14, 16, 19, 21, 23, 25-31, 14, 20, 32, 43, 52, 
76, 36-41, 44-49, 50^ 2> 51, 53i 2, 67, 68, 75, 77, 113, 114 
(omitting 83, 109); N = 46. 



Y = 252.96 + 35.175 



X 



R = 0.46873 



Significance = 0.0010 



VARIABLE 



MINIMUM 



MAXIMUM 



MEAN 



STANDARD 
DEVIATION 



Depth [meters] 
% Organic 

Mixed Uptake [yg02/gia sed/5 min] 
Average Core Uptake [yg02/cm2/hr] 
[mixed uptake] 

Ratio 

[100(% organic)] 



7. 


24. 


15.3 


5.0 


4.22 


10.19 


7.735 


1.558 


271. 


754. 


525.0 


116.9 


4.99 


8.29 


6.635 


0.767 



0.42 



1.03 



0.692 



0.144 



Regions B and Bl. Stations 15, 17, 22, 33-35, 78, 79, 84, 85, 101-103 
(omitting 81, 104, 108); N = 14. 



Y = 54.403 + 23.856X 



R = 0.53262 



Significance = 0.0499 



VARIABLE 



MINIMUM 



MAXIMUM 



MEAN 



STANDARD 
DEVIATION 



Depth [meters] 
% Organic 

Mixed Uptake [yg02/gni sed/5 min] 
Average Core Uptake [yg02/cm2/hr] 
[mixed uptake] 

Ratio 

[100(^ organic)] 



9. 


23. 


16.9 


5.0 


1.28 


8.76 


5.208 


2.09 


48. 


315. 


178.6 


93.6 


5.52 


8.33 


6.780 


0.801 



0.10 



0.53 



0.350 



0.120 
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TABLE 5. (Continued). 



Regions CI, C2, and C3. Stations 492, 53^, 54, 55]l 2"58, 63, 64, 69, 74, 80]_ 2> 
82i 2> 832, 91, 92, 94, 96, 110, 112 (omitting 14, 72, 81, 83, 153); N = 23. 



Y = -475.85 + 95..71X 



R = 0.57601 



Significance = 0.0040 



VARIABLE 



MINIMUM 



MAXIMUM 



MEAN 



STANDARD 
DEVIATION 



Depth [meters] 
% Organic 

Mixed Uptake []ig02/gni sed/5 min] 
Average Core Uptake [yg02/cin^/hr] 
[mixed uptake] 

Ratio 

[100(% organic)] 



17. 


24. 


21.0 


1.9 


7.37 


10.76 


9.388 


0.936 


171. 


706. 


422.7 


155.5 


3.20 


5.74 


4.531 


0.769 



0.20 



0.75 



0.447 



0.144 



Region D. Stations 86, 87, 97, 111, 115, 116, 118, 120, 122, 128, 129, 131, 
136, 138, 143, 145, 147-149, 152, 154, 155 (omitting 20, 32, 43, 72, 99, 100, 
119, 153, C-9); N = 22. 



Y = -79.189 + 27.178X 



R = 0.77617 



Significance 0.0000 



VARIABLE 



MINIMUld 



MAXIMUM 



MEAN 



STANDARD 
DEVIATION 



Depth [meters] 
% Organic 

Mixed Uptake [yg02/gm sed/5 min] 
Average Core Uptake [ijg02/cm2/hr] 
[mixed uptake] 

Ratio 

[100(% organic)] 



9. 


29. 


20.0 


5.5 


1.12 


7.52 


3.536 


1.778 


15. 


218. 


95.3 


62.3 


1.68 


4.54 


3.170 


0.755 



0.05 



0.53 



0.270 



0.117 



Region Eo Stations 123, 124, 127, 132, 133, 135, 137, 139-142, C-5, C-7 
(omitting 134, 150); N = 13. 

Regression not meaningful. 



VARIABLE 



MINIMU>{[ 



MAXIMUM 



MEAN 



STANDARD 
DEVIATION 



Depth [meters] 
% Organic 

Mixed Uptake [yg02/gfli sed/5 min] 
Average Core Uptake [\ig02/cxnr/hr] 
[mixed uptake] 

Ratio 

[100(% organic)] 



27. 


62. 


42.5 


9.7 


3.64 


9.32 


7.196 


1.750 


31. 


446. 


162.7 


129.1 


1.29 


4.38 


2.559 


0.970 



0.06 



0.54 



0.222 



0.161 
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TABLE 6. Oxygen uptake of mixed sediment related to percent organic matter for 
regions determined by ratio of mixed sediment oxygen uptake rate to percent 
organic matter in Lake Erie, 1963. Y = yg O2 consumed/gm sediment/ 
5 minutes; X = % organic matter. 



Region 1. Stations 5, 6, 9, 10, 12, 14, 19 (omitting 18); N = 7. 

Y = -209.12 + 125. 77x R = 0.97122 Significance 0.0003 



VARIABLE 



MINIMUM 



MAXIMUM 



MEAN 



STANDARD 
DEVIATION 



Depth [meters] 
% Organic 

Mixed Uptake [ug02/gm sed/5 min] 
Average Core Uptake [yg02/cm^/hr] 
[mixed uptake] 

Ratio 

[100(% organic)] 



9. 


13. 


10.7 


1.3 


4.22 


7.20 


5.969 


1.080 


340. 


742. 


541.6 


139.9 


4.99 


7.97 


6.350 


0.982 



0.79 



1.03 



0.897 



0.081 



Region 2. Stations 2-4, 8, 11, 13, 15, 16, 20-23, 25-32, 48, 49^, 50^ 2> 
51-55]_ 2> 72, 114 (omitting 17, 53); N = 32. 

Y = -72.798 + 79.972X R = 0.96044 Significance = 0.0000 



VARIABLE 



MINIMUM 



MAXIMUM 



MEAN 



STANDARD 
DEVIATION 



Depth [meters] 
% Organic 

Mixed Uptake []jg02/gni sed/5 min] 
Average Core Uptake [yg02/cm2/hr] 
[mixed uptake] 

Ratio 

[100(% organic)] 



7. 


22. 


13.9 


4.0 


0.42 


10.37 


6.827 


2.573 


8. 


754. 


473.2 


214.3 


2.56 


8.29 


6.235 


1.373 



0.06 



0.85 



0.636 



0.195 
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TABL£ 6. (Continued). 



Region 3. Station 33-36, 38-41, 46, 47, 56-58, 68, 69, 1^,-lh, 78-82, 84, 85, 
110-112 (omitting 37, 43, 83, 113); N = 30. 



Y = -95.766 + 57.698X 



R = 0.90501 



Significance 0.0000 



VARIABLE 



MINIMU14 



MAXIMUM 



MEAN 



STANDARD 
DEVIATION 



Depth [meters] 
% Organic 

Mixed Uptake [yg02/gm sed/5 min] 
Average Core Uptake [iig02/cm2/hr] 
[mixed uptake] 

Ratio 

[100(% organic)] 



14. 


24. 


20.5 


2.7 


2.56 


10.76 


8.177 


2.224 


77. 


648. 


376.0 


141.8 


3.43 


8.33 


5.860 


1.361 



0.25 



0.60 



0.448 



0.085 



Region 4. Stations 44, 45, 66, 77; N = 4. 
Y = -6.8055 + 68.745X R = 0.97991 



Significance 0.0201 



VARIABLE 



MINIMULI 



MAXIMUM 



MEAN 



STANDARD 
DEVIATION 



Depth [meters] 
% Organic 

Mixed Uptake [lJg02/gm sed/5 min] 
Average Core Uptake [yg02/cm^/hr] 
[mixed uptake] 

Ratio 

[100(% organic)] 



22. 


24. 


22.5 


1.0 


7.37 


9.02 


8.365 


0.708 


505. 


623. 


568.3 


49.7 


5.98 


6.33 


6.115 


0.157 



0.66 



0.69 



0.680 



0.014 



Region 5. Stations 78, 79, 81, 84-87, 99-104, 107, 108, 110, 111, 114-118, 
120, 122, 124, 127-129, 131, 132, 136, 138, 143, 147-153 (omitting 109); N = 39. 



Y = -24.174 + 41.064X 



R = 0.88539 



Significance = 0.0000 



VARIABLE 



MINIMU^L 



MAXIMUM 



MEAN 



STANDARD 
DEVIATION 



Depth [meters] 
% Organic 

Mixed Uptake []-ig02/gm sed/5 min] 
Average Core Uptake []jg02/cm^/hr] 
[mixed uptake] 

Ratio 

[100(% organic)] 



9. 


46. 


21.2 


8.1 


0.25 


7.37 


3.512 


1.881 


12. 


313. 


120.0 


87.3 


0.82 


8.33 


4.127 


1.806 



0.11 



0.76 



0.333 



0.134 
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TABLE 6. (Continued) 



Regions 5i and 62- Stations 63, 64, 91, 92, 94, 96, 123, 133, 135, 137, 
139-142, 145, 154, C-5 (omitting 97-134) N = 18. 

Significance = 0-0036 



Y = -199.00 + 47.561X 



R = 0.70091 



VARIABLE 



MINIMUM 



MAXIMUM 



MEAN 



STANDARD 
DEVIATION 



Depth [meters] 20. 55. 32.4 11.2 

% Organic 4.86 9.60 7.868 1.385 

Mixed Uptake [yg02/gm sed/5 min] 31. 355. 172.7 102.3 

Average Core Uptake [yg02/cm2/hr] 1.29 4.98 3.182 1.186 
[mixed uptake] 

Ratio 0.06 0.38 0.209 0.109 

[100(% organic)] 



Region 61. Stations 63, 64, 91, 92, 94, 96 (omitting 92); N = 6. 
Regression not meaningful. 



VARIABLE 



Depth [meters] 
% Organic 

Mixed Uptake [yg02/gm sed/5 min] 
Average Core Uptake [yg02/cm2/hr] 
[mixed uptake] 

Ratio 

[100(% organic)] 



MINIMUM 



MAXIMUM 



MEAN 



0.20 



0.38 



0.292 



STANDARD 
DEVIATION 



20. 


24. 


21.7 


1.5 


7.96 


9.60 


8.753 


0.637 


171. 


355. 


256.7 


67.9 


3.20 


4.98 


4.345 


0.660 



0.071 



Region 62. Stations 123, 133, 135, 137, 139-142, 145, 154, 155, C-5; N = 12. 
Regression not meaningful. 



VARIABLE 



MINIMUM 



MAXIMUM 



MEAN 



STANDARD 
DEVIATION 



Depth [meters] 
% Organic 

Mixed Uptake [iig02/gm sed/5 min] 
Average Core Uptake [yg02/cm2/hr] 
[mixed uptake] 

Ratio 

[100(% organic)] 



27. 


55. 


37.8 


10.0 


4.86 


9.32 


7.425 


1.462 


31. 


327. 


130.7 


91.2 


1.29 


3.72 


2.601 


0.933 



0.06 



0.37 



0.168 



0.103 
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DISCUSSION 

Organic Deposition 

Deposition patterns of organic matter depend primarily on source locations 
and circulation. Due to the highly productive nature of Lake Erie, the water 
column is a very important source of organics. Productivity increases in a 
westerly direction. This phenomenon is attributed to several factors, 
including a high degree of mixing and nutrient regeneration in the shallow 
western basin. Additionally, there is a large influx of allochthonous 
nutrients from major rivers, especially the Detroit River, discharging into 
western Lake Erie. 

Keeping these sources in mind, deposition will also depend on how far and 
in which direction this organic matter is transported. In the western basin 
and along the southern shore, a generally eastward flow predominates (Hartley 
1968), with a strong southerly current just east of the island area (Figure 9). 

In the deeper portions of the lake, a distinct two-layer circulation 
pattern exists, often in opposing directions (Hartley 1968). This is 
reasonable when the massive eastward surface flow (Figure 10), resulting from 
the prevailing westerly winds, is considered to require a return flow upwind. 
Seiches thus produced could be an important factor in resuspending sediments 
below effective wave depths. 

Assuming bottom currents as having the predominant effect on deposition 
and resuspension, it can be seen that % organic matter (Figure 3) corresponds 
well with Hartley's circulation pattern (Figure 9). Areas of high current 
intensity show little deposition (e.g. Pelee Passage), as well as those shallow 
areas exposed to wave action. Organic sinks occur in the middle of gyres (e.g. 
south of Pte. Aux Pins), in areas of opposing currents (e.g. north of Sandusky, 
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northeast of Cleveland), or in the deep hole of the eastern basin • The broad 
tongue of high organic content extending east along the south shore of the 
central basin implies a wider eastward flow than Hartley's circulation pattern 
suggests. The large region of low organic content in the northeast corner or 
the central basin is possibly a result of low input rather than resuspension. 

Despite these relationships to wave energetics, % organic matter had poor 
correlation to depth. Indeed, the entire western basin is not in equilibrium 
with respect to wave energy and sediment particle size. The presence of these 
fine particulates indicates that input overwhelms redeposition processes 
(Thomas et al. 1976) . 

Mixed Sediment Oxygen Uptake Rates 

The mixed sediment oxygen uptake rate is considered to represent the 
maximiam possible oxygen demand exerted by the sediments on the overlying water. 
Naturally, such optimal conditions for oxidation would not exist except in 
extreme storm surges and the resulting seiches. This extremely fast (within 
5 minutes) uptake of oxygen was probably a chemical oxidation of the organic 
matter and any other reduced substances present, such as iron and manganese 
ions . 

In spite of the artificiality of this method, the strong correlation to % 
organic matter suggests its usefulness as an indicator of oxygen stress due to 
sediments. Also, the relationship between these has an intriguing stipulation. 
While high and low uptake indicate high and low organic content, respectively, 
high % organic matter does not necessarily indicate high oxygen uptake. 
Possibly, older organic matter exists that has already undergone a degree of 
oxidation and therefore exerts a much lower oxygen demand. 

With this relationship in mind, an attempt was made to separate the lake 
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into regions of similar organic type. The cluster analysis was done to provide 
a means of dividing the lake into these regions based on all variables 
measured: depth, % organic matter, mixed sediment oxygen uptake, all core 
oxygen uptake rates, and the ratio value (Figure 6). The ratio value was 
devised to get an index of the oxidizability of the organic matter at each 
station individually, as well as to allow a more refined separation into 
regions based only on those variables directly of interest (Figure 8). 

The graphs of regions thus determined show a characteristic % organic 
matter coupled with a distinctive rate of oxygen uptake. The slope of the 
regression line is of particular interest in that it assesses the relative 
oxygen demand of the organic matter in the sediments. 

Especially for the regions determined by ratio similarity (Figure 8), 
there exists an easterly decrease in regression slope through Regions 1, 2, 4, 
3, 6, and 5 (Table 6). This indicates an easterly decrease in organic 
oxidizability and supports the theory that much of the organic matter 
originates in the western basin. As it is carried eastward by the previously 
mentioned circulation, this organic matter is oxidized and exerts a lesser 
oxygen demand upon settling. 

It can be seen that in Regions 1 and 2 the organic matter present exerts a 
high degree of oxygen demand, although the regions are not always high in 
organic content. Region 3, on the other hand, is further removed from the 
source and the organic matter has a moderate degree of oxidizability. Both 
Regions 2 and 3 correspond well with circulation patterns, trailing off to the 
east on both the north and the south coastlines. Region 4 corresponds to the 
middle of the gyre proposed by Hartley (1968), and is probably a sink 
characterized by high organic content. This is adequate cause for medium high 
oxygen uptake despite the relative remoteness from sources. In both Regions 5 
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and 6, there is low oxygen uptake due to the distance from organic sources. 
There are, however, distinctions based on the % organic matter; notably, the 
eastern basin (Region 62) has a very high organic content, while Region 5 
(surrounding it) behaves similarly but has much less organic matter present. 

The cluster analysis does not yield such a clear picture. This is because 
the analysis is much more encompassing, perhaps too much so, giving 
disproportionate consideration to core uptake values. Nevertheless, the 
stations thus clustered did exhibit fairly tight clumping when mixed uptake was 
plotted against % organic matter. 

The least square regressions done on these graphs were not quite as 
significant as the ratio plots (Table 5). With the exception of the Region A 
and C reversal, however, there did exist an easterly decrease in slope, the 
significance of which was discussed previously. Also, the easterly trailing 
off of regions is again seen on the north and south coastlines (Figure 6). 
The eastern basin is equally well delineated. Regions A and C extend much 
further east than the corresponding ratio Regions 2 and 3. 

Gore Oxygen Uptake Rate 

The core uptake rate, as opposed to mixed sediment uptake, was supposed to 
represent a minimal oxygen demand the sediment could exert on the overlying 
water. In the process of centrifuging twice, the normally loose sediment was 
compacted, eliminating interstitial spaces. 

The slower uptake, compared to mixed sediment uptake of oxygen, was most 
likely dominated by a biological oxygen demand due to bacterial respiration. 

The wide range of values obtained for the 2-hour incubations compared to 
the close agreement between the 4- and 6-hour incubations suggests that the 
former had insufficient time for the uptake rate to equilibrate, while the good 
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correlation between the latter two indicates that 4 hours was probably an 
adequate incubation duration. This is supported by the observation that, while 
none of the core measurements correlated well with the other variables, the 
4-hour incubation did so better than the others. 

The high correlation coefficient between each incubation period and the 
average core oxygen uptake demonstrates that the method was consistent. 
Unfortunately, the lack of correlation to other parameters made it impossible 
to evaluate this procedure as a measurement of oxygen demand exerted by the 
sediments . 

BOD Determination 

Unfortunately, BOD determinations, as well as other physical and chemical 
parameters, could not be included in the analysis because data were not present 
for all of the stations considered. Howeve^r, based on a 5-day BOD of 
approximately 1 mg/L/5 days, it would require 46 days to deplete the 
hypolimnion of oxygen, while this actually occurred in about 20 days after 
stratification. This indicates that the BOD of the overlying water is 
insufficient to be the major factor in oxygen depletion (Carr et al. 1963). 

CONCLUSION 

The results of this research are primarily consistent with the conclusions 
of Project Hypo (Burns and Ross 1972c). Sediment uptake is seen to be more 
than capable of depleting the hypolimnion of oxygen. This is made evident by 
comparing aerial uptake rates. In situ measurements of 0.40, 0.31, and 0.39 gm 
02/m2/d£Ly were obtained by Lucas and Thomas (1971), Blanton and Winklhofer 
(1972), and Burns and Ross (1972c), respectively. The average core uptake rate 

41 



of 4.987 yg/cm^/hr corresponds to 1.20 gm 02/m^/day. This suggests that, 
rather than representing minimal uptake as previously assumed, core rates are 
in excess of those occurring naturally. Apparently the disturbance of the 
oxygen gradient over sediment depth brought reduced substances to the 
sediment-water interface and was sufficient cause for increased uptake in spite 
of compaction. 

Unfortunately, the values obtained by the mixed sediment method are not 
comparable to other studies, in that oxygen uptake is expressed per gram 
sediment as opposed to aerial or volumetric uptake. 

Menon et al . (1972) reported that bacterial oxygen uptake was only 52.6% 
of theoretical maximum. This suggests that oxygen could have been at limiting 
concentration at the sediment interface and explains the increased oxygen 
uptake observed in this study by sediments exposed to oxygen-saturated water. 

While BOD values were shown to be insufficient to account for oxygen 
depletion in the hypolimnion. Burns and Ross (1972d) more clearly illustrated 
the importance of the sediments in oxygen depletion. The estimated 328 x 10^ 
moles O2 consumed by the sediments of the central basin was 81% of the total 
oxygen consumed from the hypolimnion. Also, 96% of the aerobic heterotrophic 
bacteria in the hypolimnion was found on the bottom. Approximately 88% of the 
oxygen consumed is due to this bacterial respiration and oxidation of organic 
matter; the remaining 12% is consumed by reduced metallic species. 

The assertion that there is a large flux of organic matter from west to 
east is supported by Burns et al. (1976). Additionally, the depletion pattern 
reported by Burns and Ross (1972c) is in agreement with this study, the western 
third of the lake and the southern coastline of the central basin having the 
highest oxygen uptake rates. 

While this research did not yield oxygen uptake values comparable with 
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those measured in situ by other investigators, it did produce information that 
clarifies oxygen depletion processes in Lake Erie. The use of mixed sediment 
uptake versus % organic matter regression is instrumental in determining the 
degree to which organic matter in the sediment consumes oxygen. 
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